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ABSTRACT 

We present a simple method for tracing the spatial distribution and predicting 
the physical properties of the Warm-Hot Intergalactic Medium (WHIM) , from 
the map of galaxy light in the local universe. Under the assumption that 
biasing is local and monotonic we map the ~ 2 /i~^Mpc smoothed density 
field of galaxy light into the mass density field from which we infer the spatial 
distribution of the WHIM in the local supercluster. Taking into account the 
scatter in the WHIM density-temperature and density-metallicity relation, 
extracted from the z — Q outputs of high-resolution and large box size hydro- 
dynamical cosmological simulations, we are able to quantify the probability 
of detecting WHIM signatures in the form of absorption features in the X-ray 
spectra, along arbitrary directions in the sky. To illustrate the usefulness of 
this semi-analytical method we focus on the WHIM properties in the Virgo 
Cluster region. 

Key words: Cosmology: theory - intergalactic medium - large-scale structure 
of universe - quasars: absorption lines 



1 INTRODUCTION 

As pointed out for the first time by Fukugita, Hogan 
and Peebles (1998) the census of baryons in the low 
redshift universe is significantly below expectations. In- 
deed, the baryon fraction detected in the Lya forest at 
« ~ 2 (e.g. Ranch 1998), which agrees with the recent 
WMAP measurements of the Cosmic Microwave Back- 
ground (Bennett ei al. 2003), is more than a factor of 

2 larger than that associated with stars, galaxies and 
clusters in the universe at 2 ~ (Cen & Ostriker 1999). 

Recent studies based on large box-size, hydro- 
dynamical simulations have suggested that a significant 
fraction of the baryons at z ~ are found in a gaseous 
form, with a temperature between 10^ and 10^ K in re- 
gions of moderate overdensities 5 ~ 10 — 100 (Ostriker 
& Cen 1996; Perna & Loeb 1998; Cen & Ostriker 1999; 
Dave et al. 2001; Cen et al. 2001), hence providing a pos- 



sible explanation to the missing baryons problem. This 
Warm-Hot Intergalactic Medium (WHIM) is found to 
be in the form of a network of filaments, and could be 
observed in the spectra of background bright sources in 
the X-ray and far-UV bands in the form of absorption 
lines due to elements in high ionization state, like OVI, 
OVII and OVIII (Hellsten et al. 1998, Tripp et al. 2000, 
Cen et al. 2001, Furlanetto et al. 2004). 

The detection of WHIM, however, is challenging. To 
date the most significant detection (> 3ct) in the X-ray 
band has been made in correspondence of an absorp- 
tion line in the Chandra spectrum of the Blazar Mkn 
421 that has been interpreted as due to two interview- 
ing WHIM OVII systems at 2 = 0.011 and 2 = 0.027 
and with column densities ^ 8 x 10^* cm~'^ (Nicastro et 
al. 2004). Fang et al. (2002) have reported a detection 
of an OVIII absorption line along the sight-line toward 
PKS 2155-304 with Chandra at 2: > that, however. 
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has not been confirmed by subsequent observations per- 
formed with XMM-Newton. The identification of the ab- 
sorbers in the very local universe with WHIM structures 
in the local group/supercluster is even more controver- 
sial (Nicastro et al. 2003, Savage et al. 2003). Finally, 
emission by the WHIM could also be important and 
some detections have been claimed, either by observing 
X-ray excess in clusters of galaxies (e.g. Kaastra et al. 
2003) or by observing soft X-ray structures associated 
with galaxy overdensities (Zappacosta et al. 2004). 

Although these observations arc very encouraging 
and seem to indicate that we have just started to detect 
the WHIM, various problems are still preventing us from 
a detailed investigation of its properties. 

First of all, the detection of WHIM is particularly 
difficult and requires very sensitive UV and X-ray de- 
tectors, both for absorption and for emission processes. 
With the presently available instruments, the WHIM 
absorption lines can only be detected in very high signal- 
to-noise spectra either obtained with very long expo- 
sures of bright extragalactic sources or by observing 
variable sources like Blazars when they undergo an out- 
burst phase. In either case, expected detections will be 
so rare that a detailed study of the WHIM will have to 
wait for future observations by next generation X-ray 
satellites such as XEUS and Constellation-X (Chen et 
al. 2003, Viel et al. 2003). 

The second issue, on which we will focus in this 
work, is the difficulty in modeling the WHIM. Hydro- 
dynamical simulations of structure formation in the con- 
text of cold dark matter models have certainly provided 
a major contribution to the understanding of the physics 
of the WHIM and to the characterizing its spatial dis- 
tribution. However, in general, they are not designed to 
model the actual WHIM distribution in our universe, i.e. 
they do not tell us where wc should look for the WHIM 
in our cosmic neighborhood. The hydro-dynamical sim- 
ulations of Kravtsov et al. (2002) and, more recently 
of Yoshikawa et al. (2004) , constitute the only attempt 
to model the gas distribution in the Local Supercluster 
region to address WHIM detectability. The simulations 
by Kravtsov et al. (2002) make self-consistent use of 
constraints from the MARK HI catalog of galaxy pe- 
culiar velocity to reproduce the large-scale mass den- 
sity field within ~ 100 ft~^Mpc, including major nearby 
structures like the Local Group, the Local Supercluster, 
Virgo and Coma clusters (Klypin et al. 2003). How- 
ever, their peculiar velocity constraints are only effec- 
tive above a (Gaussian) resolution scale > 5 h~^Mpc 
which is significantly larger than the typical widths of 
the WHIM structures (Viel et al. 2003). The constraints 
of Yoshikawa et al. (2004) are derived from the redshift 
space positions of IRAS 1.2Jy galaxies and are effec- 
tive on a 5 h~^M.pc Gaussian scale. As a consequence, 
in both models the actual position of the WHIM ab- 
sorber within the resolution clement of a single numeri- 
cal simulation is essentially random. A large number of 
constrained hydro-dynamical simulations, characterized 
by independent sets of Fourier modes on sub- resolution 
scales, would then be required to produce independent 



realizations of the intergalactic gas in the local universe 
from which a probability of detecting a WHIM absorber 
at a given spatial location can be computed. 

Unfortunately, this approach is too time-consuming 
to be successfully implemented using hydro-simulations 

and to allow an accurate study of all the possible phys- 
ical processes which can influence the properties of the 
WHIM at z ~ such as the effect of different amount of 
feedback in the form of galactic winds, the distribution 
of metals and the eflect of having different temperature- 
density relations. 

For this reason we have developed an alternative 
semi-analytical modeling, still based on the same prob- 
abilistic approach, in which we use the existing results 
of hydro-dynamical simulations at 2 = to trace the 
intergalactic gas from the spatial distribution of galaxy 
light in the local supercluster. The intrinsic scatter in 
the gas vs. galaxy light relation is accounted for by 
performing independent realizations of the gas distri- 
bution in the local universe which we can produce by 
means of fast Montecarlo techniques. Another advan- 
tage of our method is that our constraints, that derive 
from mapping the galaxy light, apply down to scales 
of ~ 1 h~^Mpc, comparable to the width of WHIM 
structures and much smaller than the galaxy correlar 
tion length, hence allowing us to trace the local WHIM 
with higher resolution. Obviously, these considerations 
are valid as long eis light to g£is mapping is local and well 
modeled by our reconstruction technique, which we will 
show is indeed the case. 

The outline of the paper is as follows. In Section 2 
we will describe the galaxy catalog used in our analy- 
sis. Section 3 presents the hydro-dynamical cosmologi- 
cal simulations used to model the WHIM structures. In 
Section 4 we will describe the method to produce X-ray 
absorption spectra along an arbitrary direction in the 
sky starting from the galaxy luminosity map in the lo- 
cal universe and from the physical properties of WHIM 
structures as extracted from the hydro-dynamical sim- 
ulations. Section 5 contains a testing of the modeling, 
while an application to the local universe is presented 
in Section 6. Section 7 contains our conclusion and the 
future developments of the method proposed. 



2 THE NEARBY GALAXY LUMINOSITY 
FIELD 

The galaxy sample considered in this work has been 
extracted from the latest version of the Nearby Galax- 
ies Catalog (Tully 1988). The extraction is a cube 
with cardinal axes in supergalactic coordinates that ex- 
tend to ±1500 km s~^, centered on our Galaxy. This 
boundary fully encompasses the Virgo, Ursa Major, 
Coma I, and Fornax clusters. The database strives for 
completeness among all galaxy typos known. The all- 
sky optical surveys accumulated in ZCAT (http://cfa- 
www.harvard.edu/~huchra/zcat / zcom.htm) provides 
completeness at all but low galactic latitudes for high 
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Figure 1. Left panel: Gas density as a function of mass density (in units of tlie mean density) as extracted from the z = 
output of a ACDM simulation with a box size of 25 h~^yipc (our fiducial simulation, see Section 3). The line represents the 
median value in each bin and error bars are the 1 a values. Filled contour plots represent the number density of points in the 
simulation (number density increases by an order of magnitude at each contour level). The empty contours show the number 
density of those pixels with T > 10^ K (number density increases by an order of magnitude at each contour level). Central 
panel: Gas temperature as a function of mass density. Right panel: Gas metallicity (in units of solar metallicity) as a function 
of mass density. 



surface brightness galaxies. Though dated, the HI sur- 
vey of Fisher and TuUy (1981) provides the most com- 
plete all-sky coverage of low surface brightness irregular 
galaxies. The catalog lists angular positions, redshifts, 
and B-band luminosities. There has been compensa- 
tion for the small incompleteness associated with the 
zone of avoidance through the addition of fake sources 
through reflection of real sources at higher latitudes. In 
total, the catalog contains 1968 real sources and 106 
fake sources within the cube of diameter 3000 km s~^. 
There is correction for missing information with dis- 
tance but it is very small. The information to be used is 
the luminosity density at B-band (not the galaxy num- 
ber density). The sample has a faint luminosity clip at 
Mb = —16 -I- Slog Hq/100. The catalog is essentially 
complete for galaxies of this luminosity at the bound- 
aries of the cube along the cardinal axes. Corrections 
for lost light amount to 30% at the far corners of the 
cube. 

The analysis to be discussed is performed in real 
rather than redshift space. The finger of god veloci- 
ties in clusters have been collapsed so the clusters have 
depths comparable to their projected dimensions. Then 
distances have been assigned to groups and individual 
galaxies in accordance with the output of a Least Action 
Model (Peebles 1989; Shaya, Peebles, and TuUy 1995) 
constrained by 900 individual distance measurements. 
This model provides a reasonable description of galaxy 
flows in the neighborhood of the Virgo Cluster. 



3 HYDRO-DYNAMICAL SIMULATIONS 

The fiducial hydro-dynamical simulation used here is 
a ACDM model with the following parameters; Hq — 



WOh km s-^Mpc"^ with h = 0.67, f7om = 0.30, 
Qob ~ 0.035, flA ~ 0.70, (Jg = 0.90, and the spectral 
index of the primordial power spectrum n = 1. The 
box-size of the simulation is 25 h~^Mpc comoving on 
a uniform mesh with 768"^ cells. The comoving cell size 
is 32.6 kpc and the mass of each dark matter par- 
ticle is ~ 2 X 10^ M0 (further details can be found in 
Cen et al. 2003). The simulation includes galaxy and 
star formation, energy feedback from supernova explo- 
sions, ionization radiation from massive stars and metal 
recycling due to SNe/galactic winds. Metals are ejected 
into the local gas cells where stellar particles are located 
using a yield Y — 0.02 and are followed as a separate 
variable adopting the standard solar composition. This 
simulation can be considered a good compromise be- 
tween box size and resolution to investigate the WHIM 
properties (see discussion in Cen et al. 2001). 

The X-ray and UV background of the simulation 
at 2: = 0, which are needed to properly simulate X- 
ray absorptions, are computed self-consistently from 
the sources in the simulation box. We have checked 
that there are very small differences between the hydro- 
simulation background and the estimate from ShuU et 
al. (1999): luv = /gv(£^/13.6eV)-' ■^ with I^y = 2.3x 
10"^'^ erg cm~^ Hz~^ sr~^ s~^ for the UV background 
(which includes a contribution from AGNs and starburst 
galaxies), and Ix = Ix{E / Ex)"^-'^^ {-E / Ex), 
with I°x = 1.75 X 10"^'^ erg cm"^ Hz'^ sr"^ s~^ and 
Ex = 40 keV (Boldt 1987; Fabian & Barcons 1992; 
Hellsten et al. 1998; Chen et al. 2003), for the X-ray 
background. The difference between the amplitude of 
the UV and X-ray background of the hydro-simulations 
and the estimates above is less than a factor of 2, in 
the range 0.003 < E (keV) < 4. Among the two back- 
grounds, the X-ray one plays a major role in determining 
the statistics of the simulated spectra. We will assume 
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Figure 2. Left panel: Gas density as a function of mass density (in units of tlie mean density) as extracted from the z = 
output of a ACDM simulation of 11 fc~^Mpc (see Section 3 and Cen, Nagamine, Ostriker 2004 for more details). The lines 
represent the median value in each bin and error bars are the 1 cr values. The continuous lines in all the panels refer to a 
simulation without galactic winds while the dashed ones are for a simulation with strong galactic superwinds. Filled contour 
plots represent the number density of points in the simulation without galactic winds (number density increases by an order of 
magnitude at each contour level). Central panel: Gas temperature as a function of mass density. Right panel: Gas metallicity 
(in units of solar metallicity) as a function of mass density. 



that the amplitude of the ionizing backgrounds scales 
like (1 + 2)^, in the redshift range < z < 1 (this is 
in agreement with the evolution of the radiation field 
in the simulation). This assumption will not influence 
significantly our results: given the fact that we will gen- 
erate mock-spectra at z ~ 0. 

Along with the ionizing background, the other 
quantities extracted from the hydro simulations which 
are relevant for our work are the mass density, and the 
density, temperature and metallicity of the gas. All these 
quantities have been re-binned on a coarser 256^ cubic 
grid, with cells of comoving size 97.6 kpc. In the 
rest of the paper we will label as fiducial the simula- 
tion with a box of 25 comoving Mpc and 768"^ cells 
rebinned into a mesh of 256'^. In the following we will 
refer to these quantities as the unsmoothed fields and 
will indicate them with the upper index U. 

The gas density-mass density relation in the hydro- 
dynamical simulation is shown by the contour plot 
in the left panel of Figure 1, where we also re- 
port the median value in each density bin (contin- 
uous line) and the rms values. Similarly, the cen- 
tral and right panels show the gas temperature and 
gas metallicity as a function of mass density. The 
filled contour plots represent, in all panels, the num- 
ber density of mass density elements in a given 
bin of gas-density/mass-density, gas-temperature/mass- 
density and gas-metallicity/mass-density in the hydro- 
dynamical simulation. The empty contours represent 
the number density of gas elements with T > lO^K, i.e. 
the gas effectively responsible for absorption. We prefer 
to show these plots instead of the usual gas density- 
temperature and gas density-metallicity relation (as in 
Viel et al. 2003), because our modeling starts from the 



mass distribution and uses these relations to predict the 
WHIM distribution and properties. 

The relations shown in Figure 1 are in reasonable 
agreement with the recent results obtained by Yoshida 
et al. (2002 - their Figure 2) and also by Dave et 
al. (1999), who studied the gas cooling processes in 
the framework of Smoothed Particle Hydro-dynamics 
(SPH) simulations. We note that the scatter in temper- 
ature can be very high for values of log(l -|- 5^ > — 1, 
spanning more than one order of magnitude. The mass 
density-metallicity relation in the right panel is as ex- 
pected: low density regions are the most metal poor, 
while at large densities the metallicity can be solar. In 
this plot the scatter is extremely large for the intermedi- 
ate density regime, i.e. -1 < log(l + 5^) < 0.5, while it 
is significantly smaller for larger overdensities. This plot 
has to be compared with similar plots shown in Cen & 
Ostriker (1999b), but here we find that the scatter is 
somewhat larger than in previous simulations. This can 
be due to the different amount of feedback in the hydro- 
simulation. The effect of feedback on WHIM structures 
will be investigated in a future work. 

However, to have a first hint on what properties 
of the WHIM are sensitive to the physics implemented 
into the hydro-dynamical simulations we have analyzed 
the z = of two other simulation with a smaller box 
size and at higher resolution. These simulations have a 
box size of 11 /i~^Mpc comoving on a uniform mesh of 
432^ elements and with the following cosmological pa- 
rameters: ^om = 0.29, fioi = 0.047, r^A = 0.71,0-8 = 
0.85, n = l,h = 0.7, and are extensively described in 
Cen et al. (2003) and Cen, Nagamine & Ostriker (2004). 
The main difference between this second set of simu- 
lations and our fiducial one relies in the presence of 
galactic superwinds. These are generated as feedback 
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from star formation and are normalized to match Ly- 
man Break galaxies observations. In Figure 2 we plot 
the gas density, the gas temperature and the gas metal- 
licity vs. dark matter density relations as in Figure 1. 
Wc note that, while the gas density-dark matter density 
relation is very similar, the temperature plot and the 
metallicity one show some diflerences. In particular, it 
seems that at a fixed density the temperature of the gas 
is about a factor two smaller than in the larger box size 
one. This can be determined both by the different feed- 
back implementations and by the smaller box size which 
samples a volume which is certainly not large enough to 
contain massive structures. Another major difference is 
in the mctallicity-dark matter density diagram, where 
one can see that in the case of strong galactic winds, 
represented by the dashed line, the low density regions 
are significantly metal enriched. For the higher density 
regions with 10 < 5 < 100 the discrepancy is smaller 
but still it seems that in the presence of galactic winds 
their metallicity can be a factor ^ 2 higher than in the 
no winds case. This is a very similar trend to that found 
by Cen, Nagamine & Ostriker (2004) at higher redshift 
in the Lyman— a forest and shows that winds are effec- 
tive in polluting the low density intergalactic medium 
with metals. Another feature which is worth stressing 
is the fact that both the temperature-dark matter den- 
sity and the metallicity-dark matter density relations in 
the simulation with winds show larger scatter (about a 
factor two) when compared to the simulations without 
winds. The results with no winds (continuous line) are 
however in reasonable good agreement with our fiducial 
simulation ( Figure 1), which contains a different and 
probably less effective implementation of feedback by 
galactic winds. 

In this Section we have analyzed the properties of 
the WHIM as derived from hydro-dynamical simula- 
tions with different box-sizes and different implemen- 
tation of physical processes. Among all the effects we 
found that the presence of galactic winds and the box 
sizes are important in determining the properties of the 
WHIM in particular its temperature and metallicity. In 
the semi-analytic modeling we arc going to propose in 
the following Sections, which is of course less accurate 
than the hydro-dynamical simulations described here, 
an implementation of the different relations shown in 
Figures 1 and 2 will be straightforward, allowing us to 
produce many realizations of the local universe WHIM 
distribution. 



4 FROM GALAXY LIGHT TO WHIM 

In this Section we describe a method to infer the spa- 
tial distribution of the Intergalactic Medium from the 
observed galaxy distribution in the local universe. The 
method consists of four steps, described in the next sub- 
sections, through which the B-band luminosity density 
field obtained from the TuUy catalog is used to infer the 
spatial distribution of the dark matter. Then, from the 
dark matter density field we predict the spatial distri- 



bution, temperature and metallicity of the intergalactic 
gas in the local universe. 



4.1 Light to Smoothed Mass Mapping 

To map the observed light distribution into the mass 
distribution we have applied the inversion method of 
Sigad, Bramchini and Dekel (2000). As pointed out by 
Dekel and Lahav (1999) the relation between two ran- 
dom fields like (5l and 5m, can be described by the condi- 
tional probability distribution P{Sin\5h)- If this relation 
were deterministic, then the relation between the two 
fields would be completely specified by the mean bias- 
ing function 

{5m\Sl.) = j P(5m|<5L)<5md5m - (1) 

Under the hypothesis that the relation between 5l and 
(5m is deterministic and monotonic then the mean relar 
tion ((5m|<5L) can be obtained by equating the cumula- 
tive distribution of (5l, Cl(<5l), and Jm, Cm(<5m), at given 
percentiles: 

(5ml<5L) = CmMCL(<5L)] , (2) 

where Cm^is the inverse function of Cm. 

In general, however, the relation is not determinis- 
tic since the value of the 5l at a generic location might 
not be solely determined by 5m. In fact, other physical 
quantities, such as the gas temperature, or for exam- 
ple local shear, affect the process of galaxy formation 
and evolution (e.g. Blanton et al. 1999) and thus can 
introduce a scatter in the relation usually referred to as 
stochasticity. In presence of stochasticity, C,^^[Cl(<5l)] 
may still be a good approximation for (5m|5L) as long as 
the biasing is monotonic. Branchini et al. (2004) applied 
this method to the simulation of Cen et al. (2001) to test 
the validity of this approach. More precisely, they have 
estimated the mean biasing relation from the cumula- 
tive distribution function of the ~ 2 /i^^Mpc-smoothed 
density fields of galaxy light and mass in the hydrody- 
namical simulation. They found that the (<5m|5L) mea- 
sured in the simulation agrees well with its estimate, 
C^^[Ci,{5\.)\, within the errors. 

The method has been implemented as follows. First 
we have considered the luminosity overdensity field at 
the generic grid point position, i5l(x), and have assigned 
a mass overdensity (5i(x) = C-^[Ci^{6l{:x))]. The result 
is a mean mass overdensity field that we have perturbed 
by adding a Gaussian scatter that accounts for both 
the stochasticity and the uncertainties afl^ecting the de- 
termination of {Srn\5h)- A detailed study of the mean 
biasing function of the B-band galaxy light density of 
the Tully catalog will be presented by Branchini et al. 
(2004). Here, wc only stress that the resulting biasing 
function is monotonic, i.e. it satisfies the main require- 
ment of the light to mass mapping procedure, and is 
in good agreement with the mean B-band light biasing 
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function measured in the hydrodynamical simulation of 
Cen et al. (2001), once that sampling errors and cosmic 
scatter are accounted for. 

Using the mean biasing function measured by Bran- 
chini et al. (2004) and an estimate of the scatter around 
the mean obtained from the hydrodynamical simulation, 
we have applied our mapping procedure to the observed 
B-band galaxy light map and have obtained 30 differ- 
ent realizations of the mass density field in the local 
universe that will constitute one of the input of the re- 
construction procedure presented in this work. 



4.2 Smoothed Mass to Unsmoothed Mass 
Mapping 

The mass field 6^ obtained from the light distribution 
is smoothed on a Gaussian scale = 95 kms"'^ which 
is too large to investigate the properties of the X-ray 
absorbers in the local universe. 

The problem of how to infer an unsmoothed mass 
density field, 5]^, from a smoothed one, i.e. how to re- 
cover the correlation properties on sub-resolution scale 
from large scales constraints, is well known. A self- 
consistent solution is found when is Gaussian and the 
power spectrum of mass density fluctuations is known 
a priori (Hoffman & Ribak 1991). Unfortunately, in our 
case the nonlinear evolution has caused both 6^ and 
fields to deviate from Gaussianity and we cannot use 
conventional reconstruction techniques. 

Therefore, we have developed an alternative recon- 
struction strategy that only assumes a priori knowledge 
of the mass power spectrum which we set to be the same 
as in the hydro-simulation discussed in Section 3. This 
method consists of two steps: 

i) We start from any of the 30 realizations of 5^ and 
increase the power on sub-resolution scales by adding 
Fourier modes with random phases for wavenumbers 
larger than ks = '2i:/rs in such a way to follow the 
power spectrum of the hydro simulations and to preserve 
the original spectral shape and amplitude on wavenum- 
bers larger than ks. Motivated by the fact that the 
one point probability distribution function [PDF] of 5^ 
should be close to lognormal we add the power on the 
lognormal field 5^^ = ln(l + 5^)/ < ln(l + 5^) > -1 
rather than to 5]^. We label the new mass density field 
as 5^ to indicate that random power has been added 
to the smoothed field on sub-resolution scales. 

ii) At this point, we enforce the correct mass PDF 
by performing rank ordering between and the mass 
density field extracted from the hydro-simulation, 5^. 
We indicate the resulting mass density field as 5m''''''- 

The whole purpose of this SRR reconstruction pro- 
cedure is to obtain an unsmoothed density field, (5^'*'' , 
with a power spectrum that resembles that of the orig- 
inal hydro-dynamical simulations over the whole range 
of wavenumbers and possesses the same PDF of 5^- It is 
worth stressing that this procedure is not self-consistent 



and therefore we do not expect to reconstruct the orig- 
inal mass density field on a point-by-point basis. How- 
ever, as we will show in the next section, this recon- 
struction method works reasonably well in regions of 
enhanced density where the warm hot phase of the in- 
tergalactic gas typically resides. 



4.3 Mass to Intergalactic Gas Mapping 

As a further step we map the mass overdensity field, 
(5m'''' , into gas overdensity, temperature and metallicity. 
We do that by using the relations between these quan- 
tities determined from the hydro simulation described 
in Section 2 and displayed in Figure 1. The large scat- 
ter in these plots is highly non Gaussian. As a conse- 
quence we cannot assign any of the gas properties as a 
Gaussian deviate about the mean value of the gas den- 
sity, temperature and metallicity measured in each mass 
density bin. Instead, we take the number of elements in 
each bin of the contour plots shown in Figure 1 to be 
proportional to the joint probability of any two quan- 
tities, for example, ((5gas,5m)) and assign the gas prop- 
erties (5gas,Tgas and Zgas) to the generic mass element 
by Montecarlo sampling each of the joint probability 
functions P(5gas,5m), P(Tgas, <5^), and P(Zga.,5^),. 

This allows us to infer the spatial distribution of 
the Intergalactic Medium (IGM) in the local universe 
that is guaranteed to reproduce the correct one point 
PDF of the density, temperature and metallicity of the 
gas in the hydro-simulation. We stress that in the simple 
modeling presented here an implementation of different 
metallicity-temperature relations, motivated by physi- 
cal arguments, will be straightforward. 

It is worth noting that Montecarlo sampling of the 
probability distribution function is an approximate pro- 
cedure that has the advantage of being easy to imple- 
ment. However, the large scatter introduced by this pro- 
cedure, that refiect the one present in the probability 
functions, may spoil our ability to reconstruct the prop- 
erties of the X-ray absorbers. We will check below if this 
is indeed the case. 



4.4 Intergalactic Gas to X-ray Absorbers 
Mapping 

After having simulated the gas distribution, the pho- 
toionization code CLOUDY (Ferland et al. 1998) is 
used to compute the ionization states of metals, tak- 
ing into account the UV and X-ray background. The 
density of each ion is obtained through: ni{x) = 
nH{x)XiYzQZ/ZQ, where Xj is the ionization fraction 
of the ion as determined by CLOUDY and depends on 
gas temperature, gas density and ionizing background, 
uh is the density of hydrogen atoms, Z is the metallicity 
of the element and Yz^ is the solar abundance of the 
element. 

Among the heavy elements that can produce ab- 
sorption lines in the X-rays Oxygen is the most abun- 
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Figure 3. Left: The one point probability distribution function for the original lognormal density field extracted from hydro- 
dynamical simulation 5ln (continuous line); for the same field smoothed with a Top Hat filter of radius 2500 kms~^ (dashed 
line); for the reconstructed density field, 5^''^''^ (open squares) and for the reconstructed field smoothed on the same scale (filled 
triangles). Right: Power spectra of the same four density fields, characterized by the same symbols/line-styles. 




dant one and produces the strongest lines. In this 
work we only simulate absorption spectra for the two 
strongest transitions of the ions OVII {E — 0.57 keV, 
A = 21.61, / = 0.7) and OVIII (E = 0.65 keV, 
A = I8.97A, / = 0.42), with / the oscillator strength. 
As pointed out by Viel et al. (2003), in the temperature 
and density range typical of WHIM, these two ions have 
a recombination time smaller that the Hubble time and 



thus the approximation of photoionization equilibrium 
should be reasonable. 

Given the density of a given ion along the line-of- 
sight, the optical depth in redshift-space at velocity u 
(in km s~^) is 



© 0000 RAS, MNRAS 000, 000-000 



8 M. Viel et al. 



dy m {y)V[u-y- u y (j/) , b{y)] (3) 

where ao,/ is the cross-section for the resonant absorp- 
tion and depends on A/ and //, y is the real-space 
coordinate (in km s~^), V is the standard Voigt pro- 
file normalized in real-space, b = {2kBT/mic^y^^ is 
the thermal width and wc assume that vi = viqm 
is the peculiar velocity. Velocity v and redshift z are 
related through dX/X — dv/c, where A = Xi{l + z). 
For the low column-density systems considered here, 
the Voigt profile is well approximated by a Gaussian: 
V = iV^b)-'^ exp[-(M -y- v\{y)f/h% The X-ray op- 
tical depth r will be the sum of the source contirmum, 
which we assume we can determine, and that of cq. (3). 
Finally, the transmitted flux is simply = exp(— r). 



5 TESTING THE MODELLING 

In this Section we use the hydro simulation described in 
Section 3 to test the ability of our reconstruction proce- 
dure to recover the correct spatial distribution and phys- 
ical properties of the WHIM and the OVII and OVIII 
absorbers associated to it. 

The goodness of the first stop of our procedure, 
i.e. the ability of tracing the smoothed mass distribu- 
tion using galaxy light, has been extensively tested by 
Branchini et al. (2004). The authors have shown that 
the mapping procedure is robust and unbiased and that 
random errors in the light-to-mass mapping arc much 
smaller than the stochasticity in the biEising relation, at 
least for the smoothing scale and the density range in 
which we are interested here. Both stochasticity and 
mapping uncertainties constitute a source of scatter 
which we account for by considering 30 independent re- 
alizations of the mass density field, as anticipated in 
Section 4.1. 



5.1 The dark matter density field 

As we have outlined, our procedure of mapping into 
(5m enforces the correct PDF and power spectrum. How- 
ever, the spatial correlation properties of the recon- 
structed field, (5m^^, can differ significantly from those 
of the true field S]^. Here we use the hydro-dynamical 
simulation to estimate the errors introduced by our re- 
construction procedure and see whether they can affect 
our ability of reproducing the properties of the warm 
hot gas and X-ray absorbers. 

We start from the original mass field extracted from 
the hydro-dynamical simulation, 5^, from which we ob- 
tain a field, 6^, smoothed with a Top Hat kernel of 
radius rs = 200 kms~^. Then we apply our SRR recon- 
struction procedure to obtain a field, (5ni''^, which wc 
compare to 6^. Finally, we smooth with the same 
filter and compare it with S^. The rationale behind this 
second comparison is that the constraints imposed on 



sub-resolution scales can spoil the properties of the re- 
constructed mass density field on larger scales. 

The left panel of Figure 3 shows the one point PDF 
for the original 5ln field in the simulation (continuous 
line). The lognormal field obtained with our proce- 
dure (open squares) has, by construction, the same PDF 
as the original one. However, on the scale of smoothing, 
the PDF of the SRR-reconstructcd field (filled trian- 
gles) is different from that of the original, smoothed 
field (dashed line). The discrepancy is significant in low 
density regions and vanishes in dense environments with 
Si > 10. 

The power spectra of the original and SRR- 
reconstructed fields have also been computed and are 
shown in the right panel of Figure 3 for both the 
smoothed and unsmoothed cases. Also in this case, sig- 
nificant differences between the original and the recon- 
structed density fields only exits for fe > 2 hMpc~^ i.e. 
on scales smaller than a few Mpc. Fluctuations on a 
scale r mostly come from waves of wavelength 4r, this 
mean that fluctuations on a scale > 0.8 Mpc are not 
much affected with this technique. 

We stress hero that wo do not aim at reproduc- 
ing the power spectrum of the simulated flux perfectly. 
Our reconstruction technique is not designed to repro- 
duce the two-point correlation properties of the ab- 
sorbers that, instead, may be significantly affected by 
the large scatter introduced by the Montecarlo sam- 
pling described i Section 4.3. Nevertheless the compar- 
ison with hydro-dynamical simulations shows that dis- 
crepancies exists at an acceptable level, even in the 2 
points statistics. 

Based on these results, we claim that the errors 

involved in the SRR procedure will have small impact 
on WHIM structures that have a typical size of the order 
of 1 Mpc and are characterized by overdensities larger 
than 10. In these ranges both the PDF (left panel of 
Figure 3) and the power spectrum (right panel) of the 
reconstructed mass density field are in relatively good 
agreement with those of the original simulation. 

To have a qualitative view of the structures recon- 
structed with the SRR procedure we plot in Figure 4 
a slice around the most massive object found in the 
hydro-simulation. In the left panel we show the original 
mass density field 5]^ while in the right one wo show 
our reconstructed density field, 6^^. The results are 
as expected: the densest structures and large scale fear 
tures are reproduced correctly, while on scales smaller 
than that of smoothing the S^^' is less coherent and 
the position of the small structures is not preserved. 

Motivated by the quantitative and qualitative 
agreement between the densest structures in the simula- 
tion and in the reconstruction we will explore in the next 
section other statistics more closely related to WHIM 
modeling and detectability. 
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Figure 5. Filled contour plots represent the optical depth 
binned as a function of gas density (x-axis) and metallic- 
ity (bottom panels, y-axis) and temperature (top panels, y- 
axis). The results shown in the left panels have been obtained 
with our reconstruction method (SRR), while the right pan- 
els refer to gas in the hydro-dynamical simulation (HYDRO). 
Empty contour plots (black contours) represent the number 
density of pixels. The number density increases by an order 
of magnitude with each contour level. 

5.2 The gas producing the absorption and the 
simulated spectra 

The aim of this subsection is to compare the X-ray 
absorption features associated to the OVII and OVIII 
ions in the original hydro-dynamical simulation with 
those obtained from the full reconstruction procedure 
described in Section 4. In order to do that, we draw 
spectra from the simulation box from which we com- 
pute both the optical depth and transmitted flux due 
to OVII and OVIII absorbers and compare them with 
the same quantities obtained from the results of 30 in- 
dependent reconstructions. 

In Figure 5 the filled contour plots represent the 
optical depth in OVII binned as a function of the gas 
density and metallicity (bottom panels) and gas density 
and temperature (top panels). The left panels show the 
results of one of our reconstructed maps (SRR) while 
the results in the right panel refer to the original hydro- 
dynamical simulation. The empty contours (black line) 
represent the number density distribution of the pixels 
in the whole simulations (i.e. essentially the filled con- 
tours of Figure 1). We have explicitly checked that we 
obtain very similar results for OVIII absorbers. 

This figure shows that there is a reasonable agree- 
ment between the hydro-dynamical simulation and our 



reconstructions method at least for the optical depth of 
the simulated OVII and OVIII absorbers. This agree- 
ment indicates that our reconstruction method succeeds 
in placing the gas responsible for absorption in OVII 
and OVIII in the same environments as in the origi- 
nal hydro-simulation, i.e. within regions of quite high 
metallicity, Z > 0.1 ZQ, and temperature in the range 
10^ < T{K) < W''-^. 

In the left panel of Figure 6 we show the one point 
probability distribution function (PDF) of the transmit- 
ted flux in OVII and OVIII, while in the right panel we 
show the power spectrum of the transmitted flux. In 
both plots the triangles refer to the hydro-dynamical 
simulation while the continuous lines show the average 
over 30 independent reconstructions performed with our 
technique. The results, which are similar for both the 
OVII and OVIII absorbers, show that while the flux 
probability distribution functions are reproduced very 
well, as we expect given the constraints on the mass 
PDF, the power spectra of the transmitted flux show 
some disagreement. In the range 1 < fe ( /iMpc~^) < 5 
the flux power spectrum of the hydro-dynamical simu- 
lations is about ~ 60% higher than the reconstructed 
one. We note that this discrepancy is somewhat larger 
than the corresponding one for the mass density field 
shown in the right panel of Figure 3. Moreover, in that 
case the reconstructed power at similar wavenumbers 
was systematically larger, not smaller, than the true 
one. We regard this lack of power in the transmitted 
flux as determined by the large amount of scatter in- 
volved in our modelling, that might reduce the spatial 
correlation between the absorbers on scales between 1 
and 5 /iMpc~^. On the other hand, at wavenumbers 
larger than k — 0.3 /iMpc~^, i.e. on small scales, the 
power is reasonably well preserved in both the flux and 
the mass density field. 

In addition, we show the flux probability distribu- 
tion function and the flux power spectrum for a case 
in which any form of scatter has been switched off (i.e. 
instead of Montecarlo sampling the probability distri- 
bution functions as described in Section 4.3, we simply 
adopt the deterministic mean relations). We can see, as 
expected, that the flux 1-point pdf is very different espe- 
cially for large absorbers. This means that the strongest 
absorptions in the spectra can only be reproduced by 
taking into account the large scatter in the temperature- 
density and metallicity-density relations. This confirms, 
at least qualitatively, the results of Chen et al. 2003 and 
Viel et al. 2003. Switching off the scatter also affects 
the predicted fiux of the absorbers that systematically 
underestimates that measured in the hydrodynamical 
simulation by a factor ~ 2.5. 

In Figure 7 we show the cumulative number of 
absorbers in the hydro simulations, per unit redshift 
range as a function of their column density (bottom 
axis), i.e. the number of systems with column density 
larger than a given value, and as a function of the cor- 
responding equivalent width (top axis), for OVII (filled 
triangles) and OVIII (open triangles). The relation be- 
tween the absorber equivalent width and the column 
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Figure 6. Left: Probability distribution function (pdf) for the transmitted flux in OVII (flUed triangles) and in OVIII (empty 
triangles) for the hydro-simulation. The continuous and dashed lines show the probability distribution function for the recon- 
structed flux in OVII and OVIII, respectively. The dotted line refers to the probability distribution function for OVII in the 

case in which scatter in the reconstruction method has been switched ofi^. Right: The power spectrum of the transmitted flux 
in OVII and OVIII for the hydro-simulation and from our reconstruction procedure. The same symbols/line-styles are used. 



density is taken from Sarazin (1989) and represents a 
good approximation for column densities smaller than 
lO^^cm^^. The continuous and dashed lines represent 
the results obtained from our reconstruction method 
and refer to OVII and OVIII, respectively. The compar- 
ison with the hydro-dynamical simulation shows that 
our reconstruction method is capable of reproducing 
the correct number of systems with column density 
> 3x lO^'^cm"^. Discrepancies exist for small absorption 
systems that, however, could be detected neither with 
present nor with next generation instruments (Viel et 
al. 2003). 

As a final test, wo have quantified the ability of 
reproducing the correct number and equivalent width 
of the absorption lines in simulated X-ray spectra by 
means of a hits-and-misses statistics. For this purpose 
we have only considered the spatial distribution of the 
OVII absorbers in the hydro-simulation and in the 30 
independent reconstructions. We use only the OVII ab- 
sorbers for the reason that, for a fixed column density, 
the corresponding equivalent width in OVII is about a 
factor 2 higher than the OVIII one. This means that 
probably OVII absorbers can be detected in slightly 
lower density environments than the OVIII ones. 

We have drawn 1000 independent spectra in the 
original hydro-simulation box and along the correspond- 
ing line of sights in each of the 30 SRR reconstruc- 
tions. A positive coincidence is found when an ab- 
sorber in the hydro simulation is identified within a 
bin Av = 100 kms~^ (or Ad = 200 kms~^) measured 
along the spectrum, centered on each absorber identi- 
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Figure 7. The cumulative column density distribution func- 
tion for OVII and OVII (i.e. the number of systems per unit 
redshift with an equivalent width larger than a given value), 
as obtained from the hydro-dynamical simulations (trian- 
gles - HYDRO) and with our reconstruction method (lines 
- SRR). The symbols/line-styles are the same as in Figure 
6. The top axis shows the corresponding column density in 
OVII. 
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Figure 8. Constrained probability of finding an absorber 
in the hydro-dynamical simulation map given the detec- 
tion of a corresponding absorber in the reconstructed (SRR) 
maps. This ijrobability has been computed via hits-and- 
misses statistics using 2 velocity bins of 100 kms^^ (con- 
tinuous line) and of 200 kms~^ (dashed line). The dotted 
line refers to SRR simulations in which any form of scatter 
has been taJjen out from the modelling. Error bars repre- 
sent the scatter around the mean value of the probability in 
each column density bin. The random coincidence level for a 
All = 200 km/s is shown as the horizontal thick dot-dashed 
line. The thin dot-dashed line shows the relation between 
the column densities of the absorbers reconstructed with the 
technique described in Section 5, and the absorbers found in 
the hydro-dyuamical simulation (which has to be read on the 
right Y-axis), for coincident systems. 



ability to be higher that the random coincidence level 
(thick dot-dashed line). In this figure we also plot the 
probability of finding an absorber obtained after having 
switched off any form of scatter in the temperature- 
density and metallicity density relations (dotted lino). 
We can notice that the probability significantly decrease 
by a factor 1.3. 

However, the fact that our method is capable of 
reproducing absorption features at the correct location 
> 50% of the times does not necessarily imply that the 
correct column density of the absorbers is reproduced 
too. 

The thin dot-dashed in Fig. 8 shows a fit to the rela- 
tion between the column density of the OVII absorbers 
in the reconstructions (X a:xis) and the average column 
density of the corresponding absorbers in the hydro- 
simulation (Y axis on the right). The correlation shows 
that on average the column density of the corresponding 
system found in the hydro-simulations is ~ 20% larger 
that the simulated (reconstructed) one. This means that 
our reconstruction method tends to underestimate the 
actual detectability of the WHIM absorbers. 

Overall our results indicate that our reconstruction 
method is capable of reproducing the spatial distribu- 
tion and the physical properties of the gas responsible 
for the absorptions in OVII and OVIII and that our 
maps of the WHIM can be used to predict the location 
and detectability of these structures in the Local uni- 
verse with a success rate of at least ~ 50%, increasing 
with the column density of the absorber. 



TEMPERATURE AND METALLICITY 
MAPS OF THE LOCAL UNIVERSE 



fied in the reconstructed maps. This procedure allows 
us to evaluate the constrained probability of finding an 
absorption feature in the hydro-dynamical simulation 
associated to a simulated (SRR) absorber in the recon- 
structed spectrum, binned as a function of the column 
density of the latter. 

The results are displayed in Figure 8 where it is 
shown that the probability an OVII absorption lino in 
the reconstructed spectrum associated to an absorber 
identified in the hydro-simulations, ranges between 0.4 
and 0.6 for the smaller velocity bin and is weakly de- 
pendent on the column density. For the larger velocity 
bin, Au — 200 kms^^, the probability is larger and 
ranges from 0.5 to 1 with a stronger dependence on col- 
umn density. In this second case, all the SRR absorbers 
with column density larger than 10^^ cm~^ have a corre- 
sponding absorber in the hydro- dynamical simulations. 

We note that the result is both robust, since the 
probability changes by a factor ~ 1.5, for column den- 
sities ^ 10^^'^ cm~^, when using two different velocity 
bins Aw = 100 and 200 kms~^ and significant, since 
the error bars are small enough for the measured prob- 



In this Section we apply our method to investigate the 
spatial distribution and the physical properties of the 
WHIM in the real universe. For this purpose we ap- 
ply our reconstruction procedure described in Section 
4 to the TuUy catalog of nearby galaxies and produce 
30 independent maps of the intergalactic gas in the 
local universe that account for both the stochasticity 
on the biasing relation and the scatter in the relations 
shown in Fig. 1. As wo have explicitly chocked, switch- 
ing off the stochasticity of the bias, i.e. assuming that 
the bias is deterministic, does not modify our results 
in term of WHIM properties (WHIM properties are in 
very good agreement with the maps with stochasticity) . 
This moans that all the sources of scatter that we add 
into the model, after having produced the dark matter 
density field, i.e. the gas density, gas temperature, and 
gas-metallicity scatter, are dominant compared to the 
scatter induced by the stochasticity of the bias. 

We will focus on the general properties of the local 
WHIM and on our ability to map the OVII and OVIII 
absorbers in the local observed universe. In addition, we 
will present an example of the modelling applied to the 
Virgo Cluster region in Section 6.1. 
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The three plots of Figure 9, show the predicted rela- 
tion between our input galaxy luminosity field and the 
properties of our WHIM model, i.e. they illustrate to 
what extent the galaxy luminosity traces the properties 
of the intcrgalactic gas in the local universe. The results 
are very similar to those displayed in Figure 1, which, 
however, constitute an ingredient of our reconstruction 
scheme. The curves and the scatter are fairly similar in 
the two cases, showing that regions of high galactic lu- 
minosity 10 are typically associated to the WHIM 
structures, although the scatter in the temperature and 
metallicity of the gas is very large. In the same Fig- 
ure 9 we ovcrplot as a dashed line the same quantities 
as extracted from the fiducial hydro-dynamical simula- 
tions (the 25/h comoving Mpc model 768^ cells rebinned 
into a 256^). We note that there is in general a reason- 
ably good agreement between the SRR method and the 
hydro-dynamical simulations, showing once more the 
goodness of our reconstruction procedure. 

In Figure 10 we plot the (unconstrained) probabil- 
ity of finding a pixel with an OVII optical depth above 
a given value {nhreshoid) , calculated from the set of 30 
realizations of the local universe, as a function of the 
galaxy luminosity density in the TuUy catalog. Wo stress 
that in this case we do not look for coherent systems in 
the simulated maps of the universe, as we did in Sec- 
tion 5.2, but we simply check for a pixel in the simulated 
maps to have an optical depth in OVII above a given 
threshold. Indeed, if wo assume the optical depth to be 
small (r << 1) and constant within the system, then 
an approximate relation between equivalent width and 
optical depth is given by: W{ kms~^) = tAv, where 
Av is the spatial extent of the absorber expressed in 
kms^ . In the following we will not address the dis- 
tribution of equivalent widths (or column densities) in 
coherent systems. Rather, we keep Av fixed and equal to 
the resolution of the simulated spectra, i.e. we consider 
the distribution of the optical depth in each pixel. The 
rationale behind this statistics is to investigate the rela- 
tion between WHIM and overdensity in galaxies, and to 
quantify how likely absorption features can be detected 
as a function of the density of the environment. 

The left panel of Figure 10 shows the cumulative 
probability of a pixel with r > 0.1 as a function of 

the galaxy luminosity density in the reconstructed maps 
(thick, continuous line). Pixels with t > 0.1 have a prob- 
ability to be detected larger than 50% in environments 
with (5 > 10 and this probability becomes equal to 100% 
for S ~ 30. This means that in our model the very 
dense regions are always polluted with metals that pro- 
duce OVII absorptions. However, such a small optical 
depth will be difficult to observe given the uncertain- 
ties in the source continuum and the low signal-to-noise 
of the observed spectra, unless they are organized in 
a spatially coherent system along the line of sight. If 
we ask for stronger absorptions to take place one can 
see that pixels with an OVII optical depth above 0.5 
(2.5) arc associated with galaxy overdonsities of 5 > 16 
{5 > 48) with a probability larger than 50% (35%). In 
the right panel of Figure 10 we plot the probability map 
of an optical depth larger than 0.1 in a slice of thick- 



ness ~ 0.3 /i^^Mpc containing the Virgo cluster. The 
region around Virgo is visible as a prominent clump at 
{SGX, SGY) = (—3, 12), characterized by the presence 
(probability — 1) of OVII absorbers. The thin lines plot- 
ted in the left panel of Figure 10 shows the cumulative 
probability of a pixel with r > Tthreahoid in a spherical 
region of radius 2.5 /i~^Mpc centered around the Virgo 
Cluster. In this case the probabilities are larger than in 
the whole simulated local universe (up to a factor 2 for 
Tthreshoid = 0.5) bccausc tlic colicrcnt structure of the 
Virgo cluster is quite well reproduced by our simulated 
maps. 



6.1 Warm Hot Intergalactic Medium in the 
Virgo Cluster: a semi-analytical study 

In the previous section we have shown that a semi- 
quantitative study of the WHIM detectability in the 
local universe can be done with our technique. 

Motivated by the fact that WHIM structures could 
be detected with a significant probability around the 
Virgo cluster, we study in a more quantitative way, the 
properties of the WHIM and its detectability in this re- 
gion using our semi-analytical technique. In addition, 
we investigate the robustness of our predictions by ex- 
ploring the effect of galactic winds and of the scatter 
in the relations shown in Fig. 1 implemented into the 
model. 

In the left panel of Figure 11 we show the mass 
density field, obtained from the galaxy density field, 
smoothed with a Gaussian filter of radius 95 kms~^, 

assuming the mean biasing relation (i.e. stocliasticity 
off). As we have just pointed out, the Virgo Cluster re- 
gion could be very promising for detecting the WHIM 
component since the typical overdensity reconstructed 
with our technique could bo larger than 100 in a region 
of about (2Mpc//i)'^ around the cluster centre, while it 
drops significantly at larger distances. 

In the middle panel of Figure 11 we zoom on a 
smaller region around the Virgo cluster and we check 
how different implementations of our semi-analytical 
modelling will affect the WHIM detectability. All con- 
tours represented here enclose regions with probability 
> 0.5 of finding an OVII absorber (pixel) with opti- 
cal depth r > 0.5. In particular, since in Section 6, 
we have already explored the probability of detecting 
WHIM features with our fiducial model (25//i Mpc 768^ 
mesh rebinned into a 256"^ mesh - see Section 3) based 
on the relations shown in Figure 1, we want to explore 
the effect of modifying those relations on the predicted 
WHIM property in the regions around Virgo. We fo- 
cus on the effect of galactic winds, of the scatter in the 
different relations and on the shape of the metallicity- 
density relation. The iso-probability contours enclosed 
by a continuous line refer to a model in which the effect 
of strong galactic winds is taken into account, i.e. when 
we use the metallicity-density relation shown in the left 
panel of Figure2, rather than that of Figure 1 (while we 
keep the other relations to be those of Figure 1). 
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Figure 9. Left: Relation between galaxy density as derived from the galaxy catalog and the overdensity (left panel), temperature 
(middle) and metallicity (right) of the gas. The median values in each bin are shown as continuous lines while the error bars 
la rms value. The filled contour plots have the same meaning as in Figure 1. The dashed line represent the same information 
as extracted directly from our fiducial hydro-dynamical simulations (see Section 3.) 
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Figure 10. Left: Probability of finding an optical depth in OVII larger than a given value, TthreshoWi in the local universe as 
a function of the density in galaxies. The thin lines indicate the same probability as extracted in a smaller spherical region of 
radius ~ 2.5 /i~^Mpc around the center of the Virgo Cluster. Right: Contour plots of the probability of detecting an absorption 
with a r > 0.1 in a slice of ~ 0.3 h~^T>Apc centered around the Virgo Cluster. 



The shaded areas refer instead to a model in which 
any form of scatter in the different relations of Figure 
1 is switched off. Finally, the filled contours refer to a 
model in which the metallicity is assumed to be con- 
stant to a value of 0.3 Zq, which is the mean metallic- 
ity of our fiducial simulation. One can see that galactic 
winds, which are effective in polluting the intergalac- 
tic medium, produce a non negligible detectability of 
WHIM structures at smaller overdensities than in the 
other two cases. The fiducial model, i.e. that for which 
the relations of Figure 1 hold, produces results very close 



to the galactic winds case. In the model without scatter, 
instead, the distribution of the OVII absorbers closely 
traces that of the mass, as expected. In this case the re- 
gions containing OVII absorbers wit r > 0.5 are smaller 
than in the galactic wind case, meaning that the intrin- 
sic scatter in the different relations is effective in en- 
hancing the WHIM detectability as found by Chen et al. 
(2003) and Viel et al. (2003). The model with a constant 
metallicity set to a value of 0.3 Zq tends to produce de- 
tectable regions only in very dense environments. This 
means that a significant fraction of the absorption can 
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Figure 11. Left: Slice of ~ 0.3 Mpc/h around the Virgo Cluster representing the smoothed dark matter density field derived 
from the galaxy catalog of TuUy (1988). The box indicates the region around the Virgo cluster which is considered in the other 
two panels. Center: Map of the regions with probability larger than 0.5 of an absorber with an optical depth in OVII higher 
than 0.5 in the area around the Virgo Cluster. Iso-probability contours traced with a continuous line refer to the case of strong 
galactic winds, with the metallicity-density relation right panel of Fig. 2. The dashed regions refer to a model in which any form 
of scatter in the gas density, temperature and metallicity vs. dark matter density relations is switched off. The filled contours 
refer to a model in which the metallicity has been set to a constant value of 0.3 Z©, which is the mean metallicity of the 
simulated box. Right: Probability of having an absorption feature with optical depth in OVII higher than 0.5 as a function of 
galactic overdensities in the regions around the Virgo Cluster. The result of our fiducial model is represented by the dot-dashed 
line. The thin continuous line represents the case of a the model with galactic winds. The thick continuous line and the dashed 
line show the models with constant metallicity and no scatter, respectively. 



come from regions at higher metallicities than these (see 
Figures 1 and 2). 

In order to be more quantitative we plot in the right 
panel of Figure 11 the cumulative probability of a pixel 
with an OVII optical depth > 0.5, as a function of the 
galactic overdensity of the TuUy catalog for the vari- 
ous models explored. The probability associated to the 
model with galactic winds (thin continuous line) is sig- 
nificantly larger than for the two other models, even 
at moderate overdensities of 10. We note however that 
there is very little difference between this line and the 
dot-dashed one, which refers to our fiducial simulation, 
i.e. the one in which we have implemented the relations 
shown in Figure 1. It seems that winds are effective 
in transporting metals into the low density regions but 
they produce little effects on the denser structures in 
which WHIM absorption takes place. This means that 
our prediction for the WHIM distribution are robust 
since they are weakly affected by the presence of strong 
galactic winds. This issue will be investigated in more 
detail in a future work. In the rather unphysical case of 
pure deterministic relations between the IGM proper- 
ties, i.e. the case of no scatter (dashed line), the proba- 
bility decreases by a factor of 2 in the range of overden- 
sities 20-40. The lack of any correlation between density 
and metallicity (thick continuous line) has an even more 
dramatic effect since in this case one could hope to de- 
tect OVII absorptions features only in the few higher 
density regions (although this of course depends on the 
value chosen for the constant metallicity level). 

We stress here that the modelling presented in this 
paper offers not only predictions for observers on where 



to look at in order to detect WHIM structures but allows 
to put constraints on the state of the IGM even if OVII 
and OVIII lines are not detected along a particular line- 
of-sight. 



7 DISCUSSION AND CONCLUSIONS 

Motivated by the recent interest in the warm-hot phase 
of the intergalactic medium and by the present diffi- 
culty of its detection, we have developed a semi-analytic 
model capable of predicting the spatial distribution of 
the WHIM in the local supercluster, its physical prop- 
erties and the probability of detecting it through ab- 
sorption lines in the X-ray spectra of bright background 
extragalactic sources. 

Unlike the models previously pro- 
posed by Kravtsov, Klypin & Hoffman (2002) and by 
Yoshikawa et al. (2004) our technique does not require 
to perform a full hydro-dynamical simulation. Instead, 
since it uses the already existing outputs of numerical 
simulations, it requires much less CPU time to model 
the WHIM in the local universe from the observed high 
resolution map of galaxy light density. As a consequence 
we are able to i) trace the distribution of mass and gas 
with a resolution of ~ 1 h~^}Apc (Gaussian) and ii^ per- 
form several independent realizations to account for the 
model uncertainties, through which we can quantify the 
probability of detecting an absorption feature at a given 
redshift in an X-ray spectrum drawn along an arbitrary 
direction in the sky. 
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The most serious drawback of our technique is the 
lack of self-consistency in the reconstruction procedure 
which results in a lack of spatial coherence on sub-Mpc 
scales that prevents us from reconstructing the proper- 
ties of the IGM in the local universe on a point-by-point 
basis. However, as we have explicitly checked through 
extensive numerical tests, our model correctly predicts 
the spatial location and the physical conditions of the 
WHIM, that typically resides in regions of enhanced 
densities. We have focused on the OVII and OVHI ions 
that are responsible for the strongest WHIM absorption 
features in the X-ray spectra and have shown that our 
model correctly reproduces the statistical properties of 
the WHIM absorbers. In particular, we are able to re- 
produce the correct number of absorption systems per 
unit redshift as a function of the column density, espe- 
cially for absorbers with large equivalent width. More- 
over, we are able to predict the correct spatial location 
of the OVII absorbers 40-70% of the times, significantly 
above the random level, and up to 100% for systems 
with column density of ~ lO^^cm"^, although we some- 
what underestimate their actual equivalent width. Fi- 
nally, our semi-analytic technique is capable of repro- 
ducing the global physical properties of the IGM, quan- 
tified in terms of gas density, temperature and metallic- 
ity vs. galaxy light density relations. 

It is worth pointing out that several alternative 
routes to our reconstruction procedure are also possi- 
ble. A very interesting one, proposed by the anonymous 
referee, consists of mapping galaxy light density to OVII 
and OVHI density without any reference to the matter 
density field. The idea is to determine the distribution 
of OVII and OVHI in the hydrodynamical simulation 
and use their 1-point PDF in the inversion procedure 
described in Section 4.1 to map galaxy light to OVII 
and OVHI density in the local universe, hence avoiding 
mapping mass to intergalactic gas (Section 4.3) which 
represents the main source of scatter in our method. 
Unfortunately, as we have verified, the mean OVII and 
OVHI vs galaxy density relation flattens out and be- 
comes non monotonic when the density increases. As a 
result, it is not possible to trace OVII and OVHI in re- 
gions of enhanced densities where the stronger absorbers 
typically resides unless the inversion method of Section 
4.1 could bo suitably modified. 

Encouraged by our ability of mapping the distri- 
butions of OVII and OVHI in our local universe, we 
have focused on the region around the Virgo Cluster. We 
found a large (> 0.5) probability of an OVII absorber 
with opacity (r > 0.5) in each resolution clement con- 
tained in a quasi-spherical region of radius 2.5 h~^Mpc 
centered on the Virgo Cluster that therefore represents 
the most promising site for detecting WHIM in the lo- 
cal universe (hence confirming, at least qualitatively, the 
results of Kravtsov, Klypin & Hoffman (2002)). This 
prediction is robust, since the size of this region and 
the probability detection level do not change signifi- 
cantly when one allows for uncertainties in modeling the 
IGM properties. In addition we have studied the effect 
of strong galactic winds that are effective in polluting 
low density regions, but seems not to appreciably affect 



regions of enhanced density where OVII and OVIII ab- 
sorbers are preferentially located. We have found that 
using deterministic relations to describe the IGM prop- 
erties can be quite dangerous, as it would result in a sig- 
nificant underestimate of the occurrences of OVII and 
OVIII absorbers in low density regions (an effect that 
becomes even more dramatic when assuming an IGM 
with constant metallicity). 

This semi-analytic model is therefore able to map 
the probability of detecting absorbers associated to the 
highly ionized ions in the WHIM, and to predict the 
equivalent width of the corresponding absorption fea- 
tures in X-ray spectra. This is of considerable observa- 
tional interest since it allows us to quantify the proba- 
bility of observing absorption features in the spectra of 
all bright X-ray sources, produced by the WHIM in the 
Local Supercluster. Moreover, since our method some- 
what under predicts the column density of the absorber, 
we will also be able to provide a lower limit to the sig- 
nificance of detecting such features with a specified in- 
strument. It is worth stressing that this method offers 
the possibility to put constraints on the physical prop- 
erties of the IGM (metallicity, temperature, ionization 
state) even in the presence of non detections of OVII 
and OVIII absorbers along a particular line-of-sight and 
not only when these WHIM absorption features are de- 
tected. 

Obviously, our model is also able predict the emis- 
sion lines produced by the very ions responsible for the 

X-ray absorption features very much like in the recent 
study of Yoshikawa et al. (2004), although on a much 
more local basis. 

Our model predictions can be cross-correlated with 
existing X-ray spectra to increase the significance of 
marginal detections of absorption lines. 

It is worth stressing that while our modeling is per- 
formed in real space, absorption features are observed in 
redshift space. Redshift space distortions can be easily 
accounted for in our modeling by displacing the line- 
of-sight position of the X-ray absorbers according to 
some of the available model velocity field, like the non- 
parametric one of Branchini et al. (1999). These correc- 
tions are expected to be quite minor since the cosmic 
flow in the local supercluster is notoriously cold (Klypin 
et al. 2003). 

Our model is currently limited to a very local vol- 
ume of the universe. This is the consequence of having 
required a very dense sampling in tracing the galaxy 
light. In fact, our current implementation of the light to 
mass tracing technique requires at least one object per 
resolution element, a requirement that sets the maxi- 
mum size of the explored volume. More sophisticated 
techniques, like the one recently proposed by Szapudy 
& Pan (2004), can be successfully implemented in case 
of much sparser sampling and thus can be used to ex- 
tend our WHIM predictions beyond the boundaries of 
the local supercluster. 
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